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In this work, we have considered the flat FRW model of the universe in (n + 2)-dimensions filled 
with the dark matter (perfect fluid with negligible pressure) and the modified Chaplygin gas (MCG) 
type dark energy. We present the Hubble parameter in terms of the observable parameters Q m o, 
Q x o and Ho with the redshift z and the other parameters like A, B, C, n and a. From Stern data 
set (12 points), we have obtained the bounds of the arbitrary parameters by minimizing the \ 2 test. 
The best-fit values of the parameters are obtained by 66%, 90% and 99% confidence levels. Now 
to find the bounds of the parameters and to draw the statistical confidence contour, we first fixed 
three parameters C, n, a and then fixed the three parameters A, n, a. In the first case we find the 
bounds of (A, B) and draw the contour between them for 4D(n = 2), 5D(n = 3) and 6D(n = 4). 
In the second case we fixed three different values of A as 1, 1/3, —1/3 to find the bounds of (B, C) 
and draw the contour between them. Here the parameter n determines the higher dimensions and 
we perform comparative study between three cases : 4D (n = 2), 5D (n = 3) and 6D (n — 4) 
respectively. Next due to joint analysis with BAO observation, we have also obtained the bounds 
£^ ' of the parameters {A,B) by fixing some other parameters a and A for 4D, 5D and 6D. 

D 



9i 



I. INTRODUCTION 

In recent cosmological research work, the theoretical models and range of the cosmological parameters 
(-h ■ are tested continuously by the combination of different observational astrophysical data. Before achieving 
O ' Supernova data, it seems that the universe may be occupied by energy density which is very well distributed 
over large scales [l|. In 1992, Cosmic Background Explorer (COBE) HQ data suggested that the spectrum of 
Standard Cold Dark Matter (SCDM) should be modified and proved the necessity of existence of the cosmologi- 
cal constant. The existence of non-zero cosmological constant A (which has the equation of state w\ = —1) was 
supported in 1996 The different cosmological observation of SNela [l-Q, large scale redshift surveys fol fl~o| . 
the measurements of the cosmic microwave background (CMB) [ll], EH and WMAP [H, [l4| anticipate that our 
present universe which is expanding with acceleration, preceded by a period of deceleration. The mysterious 
observational facts were not explained by the standard big bang Cosmology with perfect fluid. Thus to integrate 
the recent prediction from observational cosmology, a modification is necessary in the matter sector of the 
Einstein Gravity The unknown candidate which is responsible for this accelerating scenario, has the property 
that the positive energy density and sufficient negative pressure, known as dark energy 

HE El- The scalar field 

or quintessence fl~7| produce sufficient negative pressure to provoke acceleration in which the kinetic term is 
dominated by the potential. From recent cosmological observations including supernova data [l8| and measure- 
ments of cosmic microwave background radiation (CMBR) [l9| it is evident that in ACDM model the Universe 
is made up of about ~ 26% matter (baryonic + dark matter) and ~ 74% of a smooth vacuum energy component. 

For z > 0.01, the TONRY data set with the 230 data points [2l[ together with the 23 points from Barris ct 
al are valid. Another data set named the "gold" sample (see H) contains 156 points, which includes the 
latest points observed by HST and this covers the redshift range 1 < z < 1.6. Recently the CMBR data (for 
recent WMAP results, see [13]) strongly support il\ + fi m = 1 for FRW universe in Einstein gravity Also in 
2007, Choudhury et al [23| showed that the best-fit value of Q. m for flat model was 0.31 ± 0.08. 

One of the most effective candidate of dark energy having positive energy density and negative pressure is 
Chaplygin gas whose EOS is given by p = —B/p |2J] with B > 0. Later, it has been generalized to the form 
p = —B/p a [25|, [26[ and thereafter modified to the form p = Ap — B/p a [27]], which is known as Modified 
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Chaplygin Gas (MCG). The Modified form of Chaplygin Gas go with the 3 year WMAP and the SDSS data 



with the choice of parameters A = 0.085 and a = 1.724 [28[ which are improved constraints than the previous 
ones (-0.35 < A < 0.025) 

The drawback of the gravitational force has been successfully explained by proposing the existence of more 
than three special dimensions [30| • Today, the existence of extra dimensions [3l[ [32| arc supported by large 
number of promising model and theories. The solutions for MCG in (n + 2)-dimcnsional FRW Cosmology are 
given in section II. Here, the table of H(z) and cr{z) is presented for different values of z. The \ 2 minimum 
test for best fit values of parameters are investigated with Stern and Stern+BAO joint data analysis in section 
III. Finally, some observational conclusions are drawn. 



II. BASIC EQUATIONS AND SOLUTIONS FOR MCG IN HIGHER DIMENSIONAL FRW 

COSMOLOGY 

We consider the Jn + 2) dimensional flat (k — 0), homogeneous and isotropic universe described by FRW 
metric is given by [33l l3~i| 

ds 2 = dt 2 - a 2 (t)[dr 2 + r 2 dx 2 n ] (1) 

where a{t) is the scale factor and 

dx 2 n = d9 2 + sirfO-LdOl + + sin 2 6 1 sin 2 6-z...sin 2 6 n -idd 2 l (2) 

The modified Einstein's field equations in higher dimension are given by 



n(n + 1) /ct N 2 



2 \a 

and 



p m + Px (3) 



a n(n — 1) /d\ 2 . .. 

a 2 \a J 

where p m is the energy density of the dark matter (prcssurclcss fluid) and p x , p x are the energy density and 
pressure of dark energy (choosing 8nG = c = 1). 

Now we consider the Universe is filled with Modified Chaplygin Gas (MCG) whose equation of state (EOS) 
is given by J27| 

Px = Ap x - — , A > 0, B > 0, < a < 1 (5) 

Px 

We also assume that the dark matter and dark energy are separately conserved. So the energy conservation 
equations in homogeneous higher dimensional cosmology are 

p m + (n + l)Hp m = (6) 

and 

Px + (n + l)H{p x + Px ) = (7) 

where H is the Hubble parameter defined as H — ^. From the first conservation equation (6) we have the 
solution of p m as 

Pm=Pm {l+zT +l (8) 
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where p m o is the present value of the density of matter and z = - — 1 is the cosmological redshift. From the 
second conservation equation (7) we have the solution of the energy density p x as 



Px 



B C(l + z \(n+l)(a+l)(A+l) 



A+l 



(A + l) 



(9) 



where C is the integration constant which can be interpreted as one of the contribution of dark energy. Now 
the equation (9) can be written as 



Px — PxO 



B 



(1 + A)C + z y n+1 )( a+ i)(A+i) 



(1 + A)C + B (1 + A)C + B 
where p X Q is the present value of the dark energy density. 



(10) 



III. OBSERVATIONAL DATA ANALYSIS MECHANISM 



We now investigate the expected bounds of the theoretical parameters by x 2 statistical best fit test with the 
basis of H(z)-z (Stern) [35| and Stern+BAO (36l - l4l| joint data analysis. We also determine the statistical con- 
fidence contours between any two parameters of of MCG in FKW higher dimensional cosmology. To investigate 
the bounds of model parameters here we consider Stern (H(z)-z) data set with 12 data of H(z)-z (Stern) given 
by H3 



z 


H(z) 


a(z) 





73 


± 8 


0.1 


69 


± 12 


0.17 


83 


± 8 


0.27 


77 


± 14 


0.4 


95 


± 17.4 


0.48 


90 


± 60 


0.88 


97 


± 40.4 


0.9 


117 


± 23 


1.3 


168 


± 17.4 


1.43 


177 


± 18.2 


1.53 


140 


± 14 


1.75 


202 


± 40.4 



Table 1: H(z) and a{z) for different values of z. 



From the solutions (8) and (10) we now express the Hubble parameter H in terms of redshift parameter z 



and two dimensionless density parameters fl 



mO 



2p„ 



n(n+l)H? s 



and fl 



x0 



2p* 



as follows: 



H{z) = H 



j'O 



B 



, (l + A)C 
(1 + A)C + B (1 + A)C + B 



(1 + *) 



(n+l)( Q +l)(A+l) 



+ Q mQ (l + z) T 



(11) 



This equation can be written in the form H(z) = HqE(z), where E(z) known as normalized Hubble 
parameter contains five model parameters A, B, C, n, a beside the redshift parameter z. Now to find the bounds 
of of the parameters and to draw the statistical confidence contour (66%, 90% and 99% confidence levels) we 
first fixed three parameters C, n, a and then fixed the three parameters A, n, a. In the first case we find the 
bounds of A, B and draw the contour between them. In the second case we fixed three different values of A as 
1, 1/3, —1/3 to find the bounds of B, C and draw the contour between them. Here the parameter n determines 
the higher dimensions and wc perform comparative study between three cases : 4D (n = 2), 5D (n = 3) and 
6D (n = 4) respectively. 
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Fig.3 



Figs. 1-3 show that the variation of B against A for a — 0.0001 and C = 1.5 in 4D, 5D and 6D respectively for 
different confidence levels. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted in 

these figures for the H(z)-z (Stern) analysis. 



A. Stern (H(z)-z) Data Analysis 

Here we analyze the model parameters using twelve data (35j | of Hubble parameter for different redshift given 
by table 1. The corresponding error a(z) also available in the table. Here we perform x 2 test as a tool to find 
the minimum values of the parameters and draw the contours with three particular confidence limit. For this 
purpose we first define the % 2 statistic with 11 degree of freedom as 

2 _ \ " (H E (Ha,A,B,C,n,a,z) - H obs ) 2 

X Stern - ^A2 [ - lZ > 

where He and H b s are theoretical and observational values of Hubble parameter at different redshifts re- 
spectively and <j{z) is the corresponding error as per Table 1. Since we are interested to determine model 
parameters, Hq is considered as a nuisance parameter and can be safely marginalized. We consider the observed 
parameters f2 m rj = 0.28, £l x o = 0.72, Ho = 72 ± 8 Kms -1 Mpc -1 and a fixed prior distribution. Here we shall 
determine the model parameters A, B, C by minimizing the \ 2 statistic. The probability distribution can be 
written as 



L = J e-^ x2s ^P(Ha)dH (13) 

where P(Hq) is the prior distribution function for Hq. As per our theoretical model of MCG the two 
parameters should satisfy the two inequalities A < 1 and B > 0. We now plot the graphs for different 
confidence levels i.e., 66%, 90% and 99% confidence levels and for three different dimensions (4D, 5D and 
6D). Now our best fit analysis with Stern observational data support the theoretical range of the parameters. 
When we fix the two parameters C = 1.5 and a = 0.0001, the 66% (solid, blue), 90% (dashed, red) and 99% 
(dashed, black) contours for (A, B) are plotted in figures 1, 2 and 3 for 4D (n = 2), 5D (n = 3) and 6D (n = 4) 
respectively and we see that A becomes negative in this case. If we fix the parameter A and a = 0.0001, the 
66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours for (B,C) are plotted in (i) figures 
4-6 for 4D, 5D and 6D respectively with A = 1, (ii) figures 7-9 for 4D, 5D and 6D respectively with A = 1/3 
and (hi) figures 10-12 for 4D, 5D and 6D respectively with A = —1/3. The best fit values of (B,C) and 
minimum values of \ 2 for different values of A = 1, 1/3, — 1/3 in different dimensions are tabulated in Table 
2. For each dimension, we compare the model parameters through the values of the parameters and by the 
statistical contours. From this comparative study, one can understand the convergence of theoretical values of 
the parameters to the values of the parameters obtained from the observational data set and how it changes 
from normal four dimension to higher dimension (6D). 
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H(z)-z data (Stern) for 6D 
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Figs. 4-12 show that the variation of B against C for a — 0.0001 in 4D, 5D and 6D respectively for different 
confidence levels for A = 1 (figs. 4-6), A = 1/3 (figs. 7-9), A = -1/3 (figs. 10-12). The 66% (solid, blue), 90% (dashed, 
red) and 99% (dashed, black) contours are plotted in these figures for the H(z)-z (Stern) analysis. 



G 



n 


A 


B 


r ■ 

Amm 


2(4D) 


-0.330 


0.419 


7.065 


3(5D) 


-0.334 


0.511 


7.154 


4(6D) 


-0.381 


0.625 


7.464 



Table 2: H(z)-z (Stern): The best fit values of A, B and the minimum values of x 2 in three different dimensions 
and for a = 0.0001 and C = 1.5. 



A 


n 


B 


C 






2(4D) 


0.694 


0.774 


11.484 


1 


3(5D) 


0.731 


0.887 


15.114 




4(6D) 


0.505 


0.754 


19.191 




2(4D) 


0.396 


0.420 


8.382 


i 

3 


3(5D) 


0.449 


0.240 


10.378 




4(6D) 


0.634 


0.466 


12.769 




2(4D) 


0.411 


0.300 


7.066 


1 
3 


3(5D) 


0.254 


0.302 


7.152 




4(6D) 


0.574 


0.668 


7.632 



Table 3: H(z)-z (Stern): The best fit values of B, C and the minimum values of \ 2 for different values of A 
in three different dimensions and for a — 0.0001. 



B. Stern + BAO Joint Data Analysis 



The Baryon Acoustic Oscillations (BAO) in the primordial baryon-photon fluid, leave a characteristic signal 
on the galaxy correlation function, a bump at a scale ~ 100 Mpc, as observed by Eisenstein et al [42j| . The peaks 
and troughs seen in the angular power spectrum arise from gravity-driven acoustic oscillations of the coupled 
photon-baryon fluid in the early Universe. The interpretation of BAO measurements are the effects of non-linear 
gravitational evolution, of scale-dependent differences between the clustering of galaxies and of dark matter and 
for spectroscopic surveys, rcdshift distortions of the clustering, which can shift the BAO features. Sloan Digital 
Sky Survey (SDSS) survey is one of the first redshift survey (46748 luminous red galaxies spectroscopic sample, 
over 3816 square-degrees of sky approximately five billion light years in diameter) by which the BAO signal 
has been directly detected at a scale ~ 100 Mpc (SDSS confirmed the WMAP results that the sound horizon 
in the today's universe). The SDSS catalog provides a picture of the distribution of matter such that one can 
search for a BAO signal by seeing if there is a larger number of galaxies separated at the sound horizon. We 
shall investigate the two parameters A and B for our model using the BAO peak joint analysis for low redshift 
(with range < z < 0.35) using standard % 2 distribution. The BAO peak parameter may be defined by 



1 T 1 dz ^ 2/3 



A -E{ Z1 )V* U Jo E(z)J (14) 
where 

n m = n m0 {i + z 1 ) 3 E{ Zl )- 2 (15) 

Here, E(z) is the normalized Hubble parameter and Z\ = 0.35 is the typical redshift of the SDSS data sample. 
This quantity can be used even for more general models which do not present a large contribution of dark energy 
at early times. Now the x 2 function for the BAO measurement can be written as in the following form 



, _ (-4-0.469) 2 
Xbao Q.017 2 ^ ' 
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Figs. 13 - 15 show that the variation of B against A for a — 0.0001 and C = 1.5 in 4D, 5D and 6D respectively for 
different confidence levels. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted in 

these figures for the Stern+BAO analysis. 



where the value of the parameter A for the fiat model (k = 0) of the FRW universe is obtained by A = 
0.469 ± 0.017 using SDSS data set [i^ from luminous red galaxies survey. Now the total joint data analysis 
(Stern+BAO) for the % 2 function defined by 

XTot = Xster7i + XbAO (1*0 

Now our best fit analysis with Stern+BAO observational data support the theoretical range of the parameters. 
In figures 13, 14 and 15, we plot the graphs of (A, B) for different confidence levels 66% (solid, blue), 90% (dashed, 
red) and 99% (dashed, black) contours for 4D, 5D and 6D respectively and by fixing the other parameters C = 1.5 
and a = 0.0001. 



n 


A 


B 




2(4D) 
3(5D) 
4(6D) 


-0.323 
-0.327 
-0.373 


0.440 
0.517 
0.632 


767.456 
767.508 
767.787 



Table 4: H{z)-z (Stern+BAO): The best fit values of A, B and the minimum values of \ 2 m three different 
dimensions and for a = 0.0001 and C = 1.5. 



IV. DISCUSSIONS 

In this work, we have considered the flat FRW model of the universe in (n + 2)-dimensions filled with the 
dark matter (perfect fluid with negligible pressure) and the modified Chaplygin gas (MCG) type dark energy. 
We present the Hubble parameter in terms of the observable parameters f2 m o, ^xo and Hq with the redshift z 
and the other parameters like A, B, C, n and a. We have chosen the observed values of f2 m o = 0.28, fl x o = 0.72 
and H = 72 Kms^ 1 Mpc^ 1 . From Stern data set (12 points), we have obtained the bounds of the arbitrary 
parameters by minimizing the % 2 test. The best-fit values of the parameters are obtained by 66%, 90% and 99% 
confidence levels. Now to find the bounds of of the parameters and to draw the statistical confidence contour, 
we first fixed three parameters C',n,a and then fixed the three parameters A,n,a. In the first case we find 
the bounds of (A, B) and draw the contour between them for 4D(n = 2), 5D(n = 3) and 6D(n = 4). In the 
second case we fixed three different values of A as 1, 1/3, —1/3 to find the bounds of (B,C) and draw the 
contour between them. Here the parameter n determines the higher dimensions and we perform comparative 
study between three cases : 4D (n = 2), 5D (n = 3) and 6D (n = 4) respectively. Finally due to joint analysis 
with Stern+BAO observational data, we find the bounds of (A, B) and draw the contour between them for 
4D(n = 2), 5D(n = 3) and 6D(?i = 4). 
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We have plotted the graphs for different confidence levels i.e., 66%, 90% and 99% confidence levels and for 
three different dimensions (4D, 5D and 6D). Now our best fit analysis with Stern observational data support 
the theoretical range of the parameters. When we fix the two parameters C = f.5 and a = 0.0001, the 66% 
(solid, blue), 90% (dashed, red) and 99% (dashed, black) contours for (A,B) are plotted in figures 1, 2 and 3 
for 4D (n = 2), 5D (n = 3) and 6D (n = 4) respectively. The best fit values of (A, B) and minimum values of 
X 2 in different dimensions (4D, 5D and 6D) are tabulated in Table 2 and we see that A becomes negative in 
this case. If we fix the parameter A and a — 0.0001, the 66% (solid, blue), 90% (dashed, red) and 99% (dashed, 
black) contours for (B,C) are plotted in (i) figures 4-6 for 4D, 5D and 6D respectively with A = 1, (ii) figures 
7-9 for 4D, 5D and 6D respectively with A = 1/3 and (iii) figures 10-12 for 4D, 5D and 6D respectively with 
A = —1/3. The best fit values of (B, C) and minimum values of x 2 for different values of A = 1, 1/3, — 1/3 in 
different dimensions are tabulated in Table 3. For each dimension, we compare the model parameters through 
the values of the parameters and by the statistical contours. From this comparative study, one can understand 
the convergence of theoretical values of the parameters to the values of the parameters obtained from the 
observational data set and how it changes from normal four dimension to higher dimension (6D). Next due to 
joint analysis with Stern+BAO observational data, we have also obtained the bounds of the parameters (A, B) 
by fixing some other parameters a and C for 4D, 5D and 6D. In figures 13, 14 and 15, we have plotted the 
graphs of (A, B) for different confidence levels 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) 
contours for 4D, 5D and 6D respectively and by fixing the other parameters C — 1.5 and a = 0.0001. The 
best fit values of {A,B) and minimum values of \ 2 f° r Stern and BAO data in different dimensions (4D, 5D 
and 6D) are tabulated in Table 4. In summary, the conclusion of this discussion suggests that for different 
dimension (4D, 5D and 6D) cosmological observation can put upper bounds on the magnitude of the correc- 
tion coming from quantum gravity that may be closer to the theoretical expectation than what one would expect. 
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